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C\J ' Abstract 

> : 

Vj^ , In the context of the super symmetrized seesaw mechanism embedded in the Min- 

imal Supersymmetric Standard Model (MSSM), complex neutrino Yukawa couplings 
can induce Electric Dipole Moments (EDMs) for the charged leptons, providing an 



additional route to seesaw parameters. However, the complex neutrino Yukawa matrix 
is not the only possible source of CP violation. Even in the framework of Constrained 
MSSM (CMSSM), there are additional sources, usually attributed to the phases of the 



trilinear soft supersymmetry breaking couplings and the mu-term, which contribute 
not only to the electron EDM but also to the EDMs of neutron and heavy nuclei. In 
this work, by combining bounds on various EDMs, we analyze how the sources of CP 



violation can be discriminated by the present and planned EDM experiments. 

> : 

1 Introduction 

The atmospheric and solar neutrino data [1] as well as the KamLAND [2] and K2K [3] results 
provide strong evidence for nonzero neutrino mass. On the other hand, from kinematical 
studies [4] and cosmological observations [5], the neutrinos are known to be much lighter 
than the other fermions. There are several models that generate tiny yet nonzero masses 
for neutrinos (see, e.g. [6]) among which the seesaw mechanism [7] is arguably the most 
popular one. This mechanism introduces three Standard Model (SM) singlet neutrinos with 
masses, Mjy, which lie far above the electroweak scale. It has been shown that for Mi > 10 9 
GeV, decays of the right-handed neutrinos in the early Universe can explain the baryon 
asymmetry of the universe [8]. In addition to this, lies at intermediate scales which 
are already marked by other phenomena including supersymmetry breaking scale, gauge 
coupling unification scale and the Peccei-Quinn scale. This rough convergence of scales of 
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seemingly distinct phenomena might be related to their common or correlated origin dictated 
by first principles stemming, possibly, from superstrings. For probing physics at ultra high 
energies which are obviously beyond the reach of any man-made accelerator in foreseeable 
future, it is necessary to analyze and determine the effects of right-handed neutrinos on the 
low-energy observables. 

The Minimal Supersymmetric Standard Model (MSSM), a direct supersymmetrization of 
the SM using a minimal number of extra fields, solves the gauge hierarchy problem; moreover, 
it provides a natural candidate for cold dark matter in the universe. For explaining the 
neutrino data within the seesaw scheme, the MSSM spectrum should be enlarged by right- 
handed neutrino supermultiplets. The resulting model, which we hereafter call MSSM-RN, 
is described by the superpotential 

W = Yt^pHSEiL? - Y* j e a pH"NiLj + \m i] N 1 N 3 - fie a(3 H^H^ (1) 

whose quark sector, not shown here, is the same as in the MSSM. Here a, (3 are SU(2) indices, 
i,j are generation indices, Ljp consist of lepton doublets (z^l, -^l)/?, an d E; L contain left- 
handed anti-leptons £f L . The superfields iVj contain anti right-handed neutrinos. Without 
loss of generality, one can rotate and rephase the fields to make Yukawa couplings of charged 
leptons (Y e ) as well as the mass matrix of the right-handed neutrinos (Mij) real diagonal. In 
the calculations below, we will use this basis. 

In general, the soft supersymmetry-breaking terms (the mass-squared matrices and trilin- 
ear couplings of the sfermions) can possess flavor-changing entries which facilitate a number 
of flavor-changing neutral current processes in the hadron and lepton sectors. The existing 
experimental data thus put stringent bounds on flavor- changing entries of the soft terms. For 
instance, flavor- changing entries of the soft terms in the lepton sector can result in sizeable 
/j, — > e7, r — > cy and r — > /ry. This motivates us to go to the mSUGRA [9] or constrained 
MSSM framework where soft terms of a given type unify at the scale of gauge coupling 
unification. In other words, at the GUT scale, we take 

Csoft = -m 2 (L\Li + E}E Z + N}N t + H j d H d + H\H u ) 

- ^m 1/2 (BB + WW + gg + K.c.) 

- ^e a/3 b Hf iHSH? + H.c.) - (AZe^HSEiL? - A^e^N^ + H.c.) 

- (^MiiNW + H.c). (2) 

Here A e = a Y e and A u = a Y u . The last term is the lepton number violating neutrino 
bilinear soft term which is called the neutrino B-term. 
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As first has been shown in [10], at lower scales, the Lepton Flavor Violating (LFV) Yukawa 
coupling Y v will induce LFV contributions to the soft masses of the left-handed sleptons. 
Consequently, the strong bounds on LFV rare decays can be translated into bounds on the 
seesaw parameters. In Sec. 4, we will discuss these bounds in detail. If we assume that 
the soft terms are of the form (2) * and Y v is the only source of LFV then mass-squares 
of left-handed sleptons can be considered as another source of information on the seesaw 
parameters. It is shown in Ref. [11] that, by knowing all the entries of the mass matrices 
of neutrinos and left-handed sleptons (both their norms and phases), we can extract all 
the seesaw parameters. However, such a possibility at the moment does not seem to be 
achievable. As a result, one has to resort to finding new sources of information on the 
seesaw parameters. 

In general, the neutrino Yukawa coupling, Y u , can possess CP-odd phases, and thus 
induces electric dipole moments (EDM) for charged leptons [12, 13]. It has already been 
suggested to extract seesaw parameters from the electron EDM, d e [14] . However, for deriving 
any information from d e we must be aware of other sources of CP violation that can give a 
significant contribution to d e . In the model we are using, there are three extra sources of 
CP violation in the leptonic sector: the physical phases of the \x parameter, the universal 
trilinear coupling a Q and the neutrino 5-term''. As first has been shown in [16], the phase 
of the neutrino S-term can induce a contribution to d e . In this paper, for simplicity, we 
will set B u — 0. The phases of ao and /i can result in comparable electric dipole moments 
for the electron, neutron and mercury. More precisely, they induce d e ~ {m e /m d )d d ~ 
{m e /m u )d u ~ e{m e /m d )d d ~ e(m e /m u )d u , where d u and d d respectively are the chromo 
electric dipole moments (CEDM) of up and down quarks which contribute to the EDMs of 
mercury (dn g ) and deuteron (dr>). In principle, the phases of a and ji can induce do which 
may be detectable in future searches [17]. On the other hand, as shown in the appendix, 
the quark EDMs and CEDMs induced by the phases of Y v are too small to be detectable in 
near future. Therefore, if complex Y u is the only source of CP violation, we expect d D to be 
too small to be detectable in the near future (do is measured with ionized deuteron which 
is depleted from electrons). Based on these observations we raise the following question: 
Considering the limited accuracy of the experiments, is it possible to discern the source of 

*In practice, confirmation of this assumption is not possible. However, this assumption will be refuted if 
we find that the flavor mixing in the right-handed sector is comparable to that in the left-handed sector. 

tin fact, apart from the phases of ao and /i there are two more sources of CP-violation: the phase of 
the CKM matrix and the QCD theta term. The contribution of the former to EDMs of charged leptons is 
negligible [15]. For the latter we assume that there is a mechanism like the Peccei-Quinn mechanism that 
suppresses the CP-odd topological term in the QCD Lagrangian. 
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the CP violation? The present paper addresses this very question. 

This paper is organized as follows. In Sec. 2, we show that there is a "novel" contribution 
to de which is proportional to mi/2, and it results from the renormalization group running 
of the trilinear couplings. As will be demonstrated in the text, the new contribution can 
dominate over those previously discussed in the literature. In Sec. 3, we first review the 
experimental bounds on the EDMs. We then review how observable EDMs of neutron and 
different nuclei are related to the EDMs and CEDMs of the quarks. In Sec. 4, we represent 
our numerical results and analyze the prospects of identifying the source of CP violation. 
Conclusions are given in Sec. 5. 

2 Contribution of Y v to EDMs 

In this section, we review the effects of complex Y v on the charged lepton EDMs which has 
been previously calculated in the literature. We also discuss a new effect which has been so 
far overlooked. In the end, we point out an unexpected suppression that occurs when we 
insert realistic values for the mSUGRA parameters. Throughout this section we will assume 
that complex Y v is the only source of CP-violation. 



As it is shown in [13], inserting LFV radiative corrections to An and m| in the diagram 
shown in Fig. (1), we obtain a contribution to the EDM of the corresponding charged lepton. 
By inserting one- loop lepton flavor violating corrections to An and m|, we obtain 

7(i) /x 2a ^ ^ (V la \ (V la V 2a \ a m a m\ m\ 

M 2 

■ M(Y u k rY u k ^YrrYr} • (-2™^^^, (3) 

— * 

where S is the spin of the lepton, V is the mixing matrix of the neutralinos, m a are the 
masses of the neutralinos and 

1 ( I — x\ + 2x iLogx L 1 — x 2 E + 2x E^ogx E \ 



g{x L ,x E ) = 



2{x E -x L f\ (l-x L f (l-x E ) 
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1 / 5 - Ax L - x\ + 2(1 + 2a: L )Loga: L \ 

2(x E -x L ) { (i-x L y )■ U 

The main contribution to the diagram shown in Fig. (1) comes from the momenta around 
the supersymmetry breaking scale (M susy ); as a result we have to insert the values of AAe 
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Figure 1: A contribution to the charged lepton dipole moments. 

and Am| at M susy by taking into account the effects of running of the effective operators 
from the scale that the right-handed neutrinos decouple down to M susy . It can be shown 
that the LFV corrections to the slepton masses remain unchanged between the two scales. 
However, lowering energy from the right-handed neutrino scale down to M susy , AAe changes 
significantly. Here, the main effect comes from the gauge interaction and we can practically 
neglect the effects of Yg on the running. The factor rjd e ~ 1.5 in Eq. (3) takes care of the 
running of AA^. 

Now, let us discuss the running of the relevant parameters from the GUT scale down to 
the right-handed neutrino scale. Let us take Mqut = 2 x 10 16 GeV and Mjv ~ Y 2 (H U ) 2 /m u . 
For Y v ~ 1 and m„ ~ 0.1 eV, we find ~ 10 14 GeV so we expect that the running of pa- 
rameters from the GUT scale down to the right-handed neutrino mass scale to be suppressed 
by Log(Mff / Mq UT ) / (IGti 2 ) ~ 0.1. Thus, we can practically neglect the running of the gauge 
and Yukawa couplings as well as the gaugino and right-handed neutrino masses in this range. 
But there is a subtlety to be noted here. Although the dominant terms of both AA^ and 
Am? are enhanced by a large log factor Log M ^ T , the effect in Eq. (3), which is given by 

Im(AA^Am~), contains only one factor of Log M ffi T . This is because the leading-log parts 
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Figure 2: The two-loop correction to Ag given by m^. Vertices marked with circles are 
Yukawa vertices and the rest are gauge vertices. F^. is the auxiliary field associated with Lj. 




Figure 3: The two-loop corrections to given by m^j. Vertices marked with ® and circles 
are Yukawa vertices and A-terms, respectively. The rest are gauge vertices. 
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of AAe and Am- have the same flavor structure J2k(Y^ l )*l^og(MQ UT /Ml)Y^ j , and thus, 
Im (AAe)ieading-iog{Am 2 L )i eadin g-i og = and the dominant contribution to Im[A^Am|] 



comes from Im 



which con- 



(AAe) leading-log( Am ^) finite and Im (AAe) f inite ( Am leading-log 

tain only one large log factor. If there is a two-loop contribution to the A^ term or mass 
matrix of the left-handed sleptons [(A^) 2 _/ oop or (Am 2 ^)2-ioop with two large-log fac- 
tors, \m[AA]- loop (Am]) l r_! L op )} and Im[(A^)^ op (Am|) 2 - Zo °P)] (here L — L indices denote 
leading-log contributions) can be comparable to Im (AAe)]^^ '(Am|)j~|° t ° p . Consequently, 
inserting the 2-loop correction to A^ and 1-loop correction to Am- (or vice-versa) in the di- 
agram shown in Fig (1), we get an effect comparable to (or dominant over) Eq. (3). The 
diagrams shown in Figs. (2) and (3) give the dominant two-loop corrections to (AAe) and 
(Am|), respectively. The leading- log parts of the diagrams are * 

= j£y t £ ViVTV" (Log^Jp) 2 (5) 

and 

(Am\U = 3m 1/2 a ^ ( Lo S^f ) ' " ( 6 ) 

Inserting these diagrams in the diagram shown in Fig (1) we arrive at the following result 
r (2 ) ( x w -2a3^ 2 v v /y lo \ /Vj a y 2B \ m 1/2 m fl m| m| 




M(Y v H yy, k] ixryyr} ■ M - 4) ^g^fj L °siiF^ ( 7 ) 

This effect had been overlooked in the literature. 

Finally, as discussed in [13], in large tan/5 domain the dominant contribution takes the 
following form: 

,7(3) _ 8a sr( Vl A ( Vla i V2a\ ^m ei m a UnP 4 2 2 4 m\ m\ 

E Im[(^T^^(^)^r] (Log^V Log^S, (8) 



where 



Ml / M 2 



,/ \ 1 ( l-x 2 L + 2x L Logx L 1 - rr| + 2a;gLogXg ' 



*Note that there are similar diagrams with B replacing W in the loops. The effects of the latter is less 
than 20% of the ones we are considering here. Such a precision is beyond the scope of this paper. 
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1 / 5 - Ax L - x\ + 2(1 + 2x L )hogx L 

2(x E - x L y \ (i - x L y 

5 - 4x E - x\ + 2(1 + 2x E )^ogx E 
+ (1 - x E Y 

Note that one should insert the value of fi at the supersymmetry breaking scale in Eq. (8). 

To evaluate the order of magnitude of the EDMs, at first sight it seems that we can 
simply set all the supersymmetric parameters to some common scale m susy and take the 
values of the functions / and h in Eqs. (3,7,8) to be numbers of order 1. However, this is not 
a valid simplification because the functions / and h rapidly decrease when their arguments 
fall below unity. In the mSUGRA model we expect the mass of the lightest neutralino to 
be smaller than that of sfermions. As a result, we expect h and g to be smaller than one. 
In section 4, we will see that this effect gives rise to a suppression by two to three orders of 
magnitude. 

3 Effects of the phases of fi and ao on EDMs 

In this section, we first review the current bounds on d e , d^, do, dn g and d n and the prospects 
of improving them. We then review how we can write them in terms of Im(/i) and Im(ao). 

• Electron EDM d e : The present bound on the EDM of electron is 

d e < 1.7 x 1(T 27 e cm at 95 % CL [18] (10) 

DeMille and his Yale group are running an experiment that uses the PbO molecules 
to probe d e . Within three years they can reach a sensitivity of 10~ 29 e cm [19] and 
hopefully down to a sensitivity of 10~ 31 e cm within five years. There are proposals 
[20] for probing d e down to 10 -35 e cm level. In sum there is a very good prospect of 
measuring d e in future [21]. 

• Neutron EDM d n : The present bound on d n [22] is 

d n < 6.3 x 10~ 26 e cm at 90% CL. (11) 

This bound will be improved considerably by LANSCE [23] which will be able to probe 
d n down to 4 x 10~ 28 e cm. 

• Muon EDM, d^: The present bound on d^ [18] is 

d^ < 7 x 10~ 19 e cm. (12) 
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There are proposals to measure down to 10 24 e cm [24]. Using the storage ring of 
a neutrino factory, measurement of down to 5 x 10~ 26 will become a possibility [25]. 

• Mercury EDM d Hg : The present bound on dn g is 



which can be improved by a factor of four [26] . 

• Deuteron EDM do'- The present bound on do is very weak; however, there are 
proposals [17] to probe do down to 



Different sources of CP- violation affect the EDMs listed above differently. As a result, in 
principle by combining the information on these observables, we can discriminate between 
different sources of CP-violation. However to perform such an analysis we must be able 
to express the EDMs in terms of Im[ao], Im[/z] and Im[Y^]. In the previous section, we 
reviewed the effects of complex Y u on d e . The effects of complex a and /i on d e are also 
well understood. However, writing d n , dn g and dp in terms of the sources of CP-violation is 
more complicated. To do so, we first have to express d n , dn g and do in terms of the EDMs 
and CEDMs of light quarks (namely, d u , dd, d s , d u , dd and d s ) and then calculate the quark 
EDMs and CEDMs in terms of Im[a ], Im[/x] and Im^]. The quark EDMs and CEDMs in 
terms of Im[a ] and Im[//] have already been calculated in the literature. In this paper we 
have used the results of Ref. [27]. As we discussed in the appendix, the effects of Imfl^] on 
the quark EDMs and CEDMs are negligible. Unfortunately, the first step (expressing d n , 
d Hg and d D in terms of the quark EDMs and CEDMs) is quite challenging. Let us consider 
them one by one. 

• d n (d q ,d q ): 

Despite of the rich literature on d n in terms of the quark EDMs and CEDMs, the 
results are quite model dependent. For example, the SU(3) chiral model [28] and QCD 
sum rules [29] predict different contributions from d u and dd to d n . Considering these 
discrepancies in the literature, in this paper we do not use bounds on d n in our analysis. 
As it is shown in [30], information on d n can help to refute the "cancelation" scenario. 
We will come back to this point later. 



d Hg \ < 2.1 x 10" i8 e cm 



(13) 



d D \ < (1 - 3) x 10~ 27 e cm. 



(14) 
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• dHg(d q ,d q ): 

There is an extensive literature on dn g [31]. In this paper, following Ref. [32], we will 
interpret the bound on dn g as 

\d d - d u \ < 2 x 1(T 26 cm. (15) 

As shown in the recent paper [33], the EDM of electrons in the mercury atom can give 
a non-negligible contribution to du g - As a result, improvements on the bound on dn g 
will not be very helpful for us to discriminate between different sources of CP- violation; 
i. e., du g also obtains a correction from complex Y v through d e . 

• d D (d q ,d q ): 

Searches for do can serve as an ideal probe for the existence of sources of CP- violation 
other than complex Y v because i) there is a good prospect of improving the bound on 
do [17]; ii) an ionized deuteron does not contain any electrons and hence we expect 
only a negligible and undetectable contribution from Y u to d D . 

To calculate d D in terms of quark EDMs and CEDMs, two techniques have been sug- 
gested in the literature: i) QCD sum rules [34] and ii) SU(3) chiral theory [28]. Within 
the error bars, the two models agree on the contribution from dd — d u which is the dom- 
inant one. However, the results of the two models on the sub-dominant contributions 
are not compatible. Apart from this discrepancy, there is a large uncertainty in the 
contribution of the dominant term: 

d D (d q , d q ) ~ -e(d u - d d ) 5tl l . (16) 

In this paper we take "the best fit" for our analysis. 



4 Numerical analysis 

In this section, we first describe how we produce the random seesaw parameters compatible 
with the data. We then describe the figs. (4-9) and, in the end, discuss what can be inferred 
from the future data considering different possible situations one by one. 

In figures (4-9), the dots marked with "+" represent d e resulting from complex Y u . To 
extract random Y v and M N compatible with data, we have followed the recipe described in 
[35] and solved the following two equations 

r]mYj^Y v {v 2 sin 2 (5)/2 = U ■ $ • M® iaa ■ $ • U T (17) 
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and 



a 


d* 



d 
6 
c 



(18) 



where v = 247 GeV, M is the mass matrix of the right-handed neutrinos, U is the mixing 
matrix of neutrinos with S13 = and $ is diag[l, e 1 ^ 1 , e 1 ^ 2 ] with random values of 0i and (f>2 
in the range (0, 2ir). Finally, M® %ag = diag[m 1 , Jm\ + Ara^ 



mf + Am^] where mi picks 
up random values between and 0.5 eV in a linear scale. The upper limit on m 1 is what 
has been found in [36] by taking the dark energy equation of state a free (but constant) 
parameter. In the above equation, r\ mu takes care of the running of the neutrino mass matrix 
from M to M susy . Since the deviation of r\ mv from unity is small [37], we have set r\ mv = 1. 

In order to satisfy the strong bounds on Br(/i — > e^y) [18] and Br(r — > /ry) [39], the 
matrix h, defined in Eq. (18), is taken to have this specific pattern with zero e/i and /it 
elements. Actually these branching ratios put bounds on (Am|) e/i and (Am|) MT rather than 
on h efl and h^. Notice that only the dominant term of Am| is proportional to h. There 
is also a subdominant "finite" contribution to Am- which is about 10% of the dominant 
effect and is not proportional to the matrix h [13]. As we saw in section 2, this finite part 
plays a crucial role in giving rise to EDMs because the dominant leading-log part cancels 
out. Nonetheless, for extracting the seesaw parameters, 20% accuracy is enough and we can 
neglect the subdominant part and take Am| proportional to the matrix h. In Eq. (18), a, 
b, c are real numbers which take random values between and 5. On the other hand, \d\ 
takes random values between and the upper bound from Br(r — > e^) [40]. To calculate 
the upper bound on \d\, we have used the formulae derived in Ref. [41]. The phase of d 
takes random values between and 2n. With the above bounds on the random variables, 
the Yukawa couplings can be relatively large, giving rise to 



Im 



Ml 



Y 

M? v 



few x 10 



and 



Im 



M 2 

rW M gut y yty 



few x 0.1. 



(19) 



(20) 



As we discussed in the end of Sec. 2, because of the presence of the rapidly changing 
functions g(xL,xn) and h(xL,xn) in Eqs. (3, 7, 8), the value of d e strongly depends on the 
values of the supersymmetric parameters. To perform this analysis we have taken various 
values of tan f3 and a and calculated the spectrum of the supersymmetric parameters along 
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the mi/2 — itlq strips parameterized in Ref. [42]. Notice that Ref. [42] has already removed 
the parameter range for which color or charge condensation takes place. 

In the figures, we have also drawn the present bound on d e [18] as well as the limits which 
can be probed in the future. The present bound is shown by a dashed dark blue line and lies 
several orders of magnitude above the d e from phases of Y v . After five years of data-taking, 
the Yale group can probe d e down to 10~ 31 e cm [19] which is shown with a dot-dashed cyan 
line in the figures. As it is demonstrated in the figures, only for large values of ao the effect 
of complex Y v on d e can be probed by the Yale group and for most of the parameter space 
the effect remains beyond the reach of this experiment. 

There are proposals [20] to use solid state techniques to probe d e down to 10~ 35 e cm 
(shown with dot-dashed yellow line in the figure). In this it can be deduced from 

the figure, we will have a great chance of being sensitive to the effects of the phases of Y v on 
d e . However, unfortunately, the feasibility and time scale of the solid state technique is still 
uncertain. 

Although for intermediate values of tan (3, the effect of the phases of Y v on d e is very low 
(< 1CT 30 e cm ), its effect can still be much higher than the four-loop effect on d e in the SM 
(the effect of the CP-violating phase of the CKM matrix) which is estimated to be ~ 1CT 38 
e cm [15]. 

In figs. (4-6) as well as in fig. 9, d e resulting from Im[/x] is also depicted. The red solid 
lines in these figures show d e from Im[/x] assuming that the corresponding du g saturates the 
present bound [26]. As it is well-known, there are uncertainties both in the value of m d [18] 
and in the interpretation of du g in terms of more fundamental parameters d u , da and d s . 
To draw this curve we have assumed m d = 5 MeV and d u — d d < 2 x 1CT 26 e cm . As it 
is shown in the figure this bound is weaker than even the present direct bound on d e . The 
purple dotted lines in figs. (4, 5, 6, 9), represent d e induced by values of Im[/x] that give rise 
to d u — d d = 2 x 1CT 28 cm (corresponding to d D = 1CT 27 e cm and d D = 5e(d d — d u )). Notice 
that these curves lie well below the direct bound on d e but the Yale group will be able to 
probe even smaller values of d e . Similarly in figs. (7,8), d e resulting from Im[ao] is depicted. 

The following comments are in order: 

1) The combination of the seesaw parameters that enter the formula for d e resulting from 
Im[Am|m 2 Am|] [see Eq. (8)] is 



Im 



yU 2 M GUT y m 2yU M GUT y 

u M 2 Y » m ? Y » Ln M 2 Y » 



~ m r Im 



V;Ln 2 ^W^ ( YjLn^^Y u 
M 2 \ M 2 



(21) 

where M is the mass matrix of the right-handed neutrinos. In contrast to this, the "new" 
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effect given in Eq. (7) is proportional to 



Im YjW^Y^Ln 





For the specific pattern of the h matrix shown in Eq. (18) (with zero e/x element) this effect 
is also given by 



For the values of supersymmetric parameters chosen in Fig. (4) (that is, sgn(/i)=+, 
tan/5 = 10, a = 0), the "new" effect is dominant and is —5 times the effect previously 
discussed in the literature. However, for a = 1000 GeV and 2000 GeV (Figs. 7,8 and 9) the 
dominant contribution is the one given by Eq. (8). 

2) In the figures, the bounds from d Hg and d D appear almost as horizontal lines. This 
results from the fact that for the m -m 1( / 2 strips that we analyze, m is almost proportional 
to mi/2- Using dimensional analysis we can write 



3) As discussed in Ref. [43], at two-loop level, the imaginary a can induce an imaginary 
correction to the Wino mass, giving rise to another contribution to the EDMs. In our 
analysis, we have taken this effect into account but it seems to be subdominant. 

In the following, we will discuss what can be inferred about the sources of CP-violation 
from d e and do if their values (or the bounds on them) turn out to be in certain ranges. 

According to the Figs. (4-6), for a = 0, any signal found by the Yale group implies that 
there are sources of CP-violation other than the phases of the Yukawa couplings. However, 
for larger values of a , the effect of Y v on the EDMs can be observed by the Yale group within 
five years. According to Figs. (7-9), for a > 1000 GeV EDMs originating from complex Y u 
can be large enough to be observed by the Yale group. Therefore, if after five years the 
Yale group reports a null result, we can derive bounds on certain combinations of seesaw 
parameters and a . At least it will be possible to discriminate between low and high a 
values. If after five years the Yale group reports a null result, we can derive bounds on the 
seesaw parameters. However, if the Yale group finds that 10~ 31 e cm < d e < 10~ 29 e cm we 




(22) 



In other words, the two effects are proportional to each other. 




where k\ and k 2 are given by the relevant fermion masses and are independent of m^. As 
a result, for a given value of d u — da, Im[/i] (or Im[a ]) itself is proportional to my 2 so d e 
will not vary with m^. 
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will not be able to determine whether d e originates from complex Y v or from more familiar 
sources such as complex ao or fi. To be able to make such a distinction, values of do down to 

10-28 _ 10 "29 

e cm have to be probed which, at the moment, does not seem to be achievable. 
If future searches for djj find dp > 10~ 27 e cm but the Yale group finds d e < 2 x 10~ 29 
e cm (this can be tested within only 3 years of data taking by the Yale group [19]), we might 
conclude that the source of CP- violation is something other than pure Im[/i] or Im[ao]; e.g., 
QCD #-term which can give a significant contribution to do but only a negligible contribution 
to d e . Another possibility is that there is a cancelation between the contributions of Im[/i] 
and Im[a ] to d e . The information on d n would then help us to resolve this ambiguity 
provided that the theoretical uncertainties in calculation of d n as well as ^ are sufficiently 
reduced. 

On the other hand, if the Yale group detects d e > 2 x 1CT 29 e cm, we will expect that 
do > 10~ 27 e cm which will be a strong motivation for building a deuteron storage ring 
and searching for d D . If such a detector finds a null result, within this framework the 
explanation will be quite non-trivial requiring some fine-tuned cancelation between different 
contributions. 

According to these figures, in the foreseeable future, we will not be able to extract any 
information on the seesaw parameters from EDMs, because even if we develop techniques 
to probe d e as small as 1CT 35 e — cm, we will not be able to subtract (or dismiss) the effect 
coming from Im[/x] and Im[ao] unless we are able to probe d u — dd at least 5 orders of 
magnitude below its present bound which seems impractical. Remember that this is under 
the optimistic assumptions that the mass of the lightest neutrino, m b and Br(r — > ej) are 
close to their upper bounds and there is no cancelation between different contributions to 
the EDMs. 

If, in the future, we realize that mi and Br(r — > ej) are indeed close to the present upper 
bounds on them and a = (a = 1000 GeV) but find d e < 10" 35 e cm (d e < 10~ 34 e cm ), 
we will be able to draw bounds on the phases of Y v which along with the information on the 
phases of the Dirac and Majorana phases of the neutrino mass matrix and the CP-violating 
phase of the left-handed slepton mass matrix may have some implication for leptogenesis. 
This is however quite an unlikely situation. 

Let us now discuss the d^. As we saw in Sec. 2, the phases of Y v manifest themselves in 
the dfj, through Im[AA^Am|] w and Im[Am|m 2 Am|] w . If a is a real number, the matrix 
Ai remains Hermitian [13]. That is the radiative corrections due to Y u cannot induce nonzero 
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Im[AA^]jj. So, we can write 

Im[A^Am|] w = Im[(A^)^(Am|) eM ] + Im[(A^) MT (Am|) T/1 ] 

and 

Im[Am|m^Am|] w = lm[(Am%) f , e m 2 e (Am 2 L ) efM } + Im[(Am|) MT m 2 (Am~) rM ]. 

The strong bounds on Br(/x — > e7) [18] and Br(r — > /ry) [39] can be translated into bounds 
on (Am|) e(U and (Am|) T/Lt as well as the corresponding elements of AA?. As a result, in the 
framework that imaginary Y v is the only source of CP-violation, we expect 

m e (Am z L ) Te 

which will not be observable even if the muon storage ring of a nu-factory is built [25]. On 
the other hand, imaginary a and ji induce ~ c? e m M /m e and allow d e to be as large as 
the experimental upper bound on it. In this case, we may have a chance of observing d^. 
Observing d^ 3> m^d e /m e will indicate that this simplified version of the MSSM is not valid. 



5 Summary and conclusions 

In this work we have studied EDMs of particles in the context of supersymmetric seesaw 
mechanism. We have examined various contributions to electron EDM induced by the CP- 
odd phases in the neutrino Yukawa matrix. Our analysis takes into account various contri- 
butions available in the literature as well as a new one, proportional to the gaugino masses, 
which is presented in Eq. (7). 

In our discussions we have first produced random complex neutrino Yukawa couplings 
consistent with the bounds from LFV rare decays and then calculated the electron EDM 
they induce along post-WMAP m — strips for given values of tan/5 and a [42]. We 
have found that, for small values of a , the new contribution (7) can be dominant over the 
other contributions from Y v that had already been studied in the literature. 

It turns out that for a realistic mass spectrum of supersymmetric particles, there is 
an extra suppression factor of 1CT 2 — 1CT 3 with which we would not encounter if all the 
supersymmetric masses were taken to be equal to each other. In figs. 4-9, the values of d e 
corresponding to different random complex Y u textures are represented by dots. For small 
values of tan/3 (tan/5 < 10) and a (a < 1000 GeV), d e induced by Y v is beyond the reach 
of the ongoing experiments [19]. Such values of d e can however be probed by the proposed 
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solid state based experiments [20]. For larger values of tan/3 and/or ao, the Yale group may 
be able to detect the effects of complex Y v on d e . As it is demonstrated in Figs. 8 and 9, for 
tan f3 = 20 and ao = 1000 — 2000 GeV, a large fraction of parameter space yields d e detectable 
by the Yale group. However, even in this case we will not be able to extract information 
on the seesaw parameters from d e because the source of CP- violation might be a and/or 
/j, rather than Y v . If the future searches for do [34] find out that do > 10~ 27 e cm then 
we will conclude that there is a source of CP-violation other than complex Y v . However, to 
prove that the dominant contribution to d e detected by the Yale group comes from complex 
Y v - hence to be able to extract information on the seesaw parameters from it- we should 
show that do < 10~ 28 — 10~ 29 e cm which is beyond the reach of even the current proposals. 
Notice that for the purpose of discriminating between complex Y v and a //i as sources of 
CP- violation, searching for dn g is not very helpful because mercury atom contains electron 
and hence du g obtains a contribution from complex Y v . That is while ionized deuteron used 
for measuring d D does not contain any electron and the contribution of complex Y u to it 
is negligible. To obtain information from d n , the theoretical uncertainties first have to be 
resolved. 

In this paper, we have also shown that for the neutrino Yukawa couplings satisfying the 
current bounds from the LFV rare decays, the electric dipole moment of muon induced by 
Y u is negligible and cannot be detected in the foreseeable future. Detecting a sizeable d M will 
indicate that there are sources of CP- violation beyond the complex Y v . 

6 Appendix 

Since do is dominantly given by d d — d u , in this section, we concentrate on evaluating d q . 
One can repeat a similar discussion for d q . 

In Sec. 2, we saw that integrating out N i: the effects of CP-violating phases appear in the 
left-right mixing of sleptons which can be evaluated to m£m susy F(Y u , Log[M^ UT / M 2 ]) / (16-7T 2 ) 2 
or for large tan/3, to m s m susy tan I3{m 2 tan 2 [3/v 2 )F'(Y u , Log[Mg [/T /M 2 ])/(167r 2 ) 2 . For ran- 
dom Y u consistent with observed m u and bounds on the branching ratios of LFV rare decay, 
the functions F and F' take values smaller than 0.1. Since quarks do not directly couple 
to the leptonic sector, the CP-violation in the leptonic sector should be transferred to the 
quark sector through one-loop (or higher loop) effective operators made of Higgs and gauge 
bosons or their superpartners. To construct such an effective operator one more factor of 
Yg is needed to compensate for the left-right mixing mentioned above. Considering the fact 
that Y T ^> Y M , Y e , the main contribution to the effective operator comes from the diagrams 
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with r and f propagating in them. So the CP-odd effective potential will be given by 



2\3 



or for large tan/5, 



m" m 2 T t a nf3F(Y„Mg[M 2 GUT /M 2 })/(lG7r 



m 2 T tan 2 [3(m 2 T tan 2 [3/v 2 )F'(Y u , Log[M GUT / M 2 }) / (16% 2 ) 3 . 



susy 

In these formula, the power of m susy , n, is determined by the dimension of the specific 
operator under consideration. 

To evaluate d g , we have to insert the CP-odd effective operator in another one-loop 
diagram. Since CEDMs mix left- and right-handed, the latter diagram should involve a 
factor of Y q . So, we can write 

7 Y q g 2 g s m 2 T t & n (3 F(Y V , Log [M 2 GUT /M 2 ]) 

and for large tan f3, 

Y q g 2 g s m 2 T tan 2 (3{m 2 T tan 2 (3/v 2 )F\Y v , Log[M GUT / M 2 ]) 



d 



167T 2 (167T 2 ) 3 m| 



susy 



As a result, we expect d d < 10~ 30 cm which cannot be observed even if the recent proposal 
[17] is implemented. We expect d u to be even smaller because Y u /Y d = m u /m d cot/? <C I. 
Notice that although d q is suppressed by a factor of m 2 tan 2 (3 / (lQir 2 m 2 susy ) , ed d can be 
comparable to d e . This originates from two facts: Y d /Y e ~ 10 and in the case of d e , as 
we discussed in Sec. 2, there is an extra suppression given by the functions g(xi,XE) and 
h(xL, xe)- If we do precise two-loop calculation of d d for a realistic SUSY spectrum, we may 
encounter similar suppression. As the above analysis show we do not expect an observable 
effect due to Y v in future searches for do and dug so it seems there is not a strong motivation 
for performing such a complicated two-loop calculation. 
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Figure 4: Electric dipole moment of the electron for a = 0, tan/3 = 10 and sgn(/x) = +. To 
draw the red solid and purple dotted lines, we have assumed that Im[/x] is the only source of 
CP- violation and have taken da — d u equal to 2 x 10 -26 cm and 2 x 10~ 28 cm, respectively to 
derive Im[/i]. To produce the dots, we have assumed that complex Y v is the only source of 
CP-violation and have randomly produced Y v compatible with the data. The blue dashed 
line is the present bound on d e [18] and dot-dashed lines show the values of d e that can be 
probed in the future [19, 20] 
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Figure 5: The same as fig. 4 for ao = 0, tan/5 = 20 and sgn(/i) = +. 
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Figure 6: The same as fig. 4 for a = 0, tan/5 = 20 and sgn(/x) = — . 
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Figure 7: Electric dipole moment of the electron for = 1000 GeV, tan/3 = 10 and 
sgn(/i) = +. To draw the red solid and purple dotted lines, we have assumed that Im[a ] is 
the only source of CP- violation and have taken d^ — du equal to 2 x 10~ 26 cm and 2 x 10~ 28 cm, 
respectively to derive Im[a ]. To produce the dots, we have assumed that complex Y v is the 
only source of CP-violation and have randomly produced Y v compatible with the data. The 
blue dashed line is the present bound on d e [18] and dot-dashed lines show the values of d e 
that can be probed in the future [19, 20] 
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Figure 8: The same as fig. 7 for a = 1000 GeV, tan/? = 20 and sgn(/x) = +. 
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Figure 9: Electric dipole moment of the electron for = 2000 GeV, tan/3 = 20 and 
sgn(/i) = +. To draw the red solid and purple dotted lines, we have assumed that lm[/j] is 
the only source of CP- violation and have taken dd — d u equal to 2 x 10~ 26 cm and 2 x 10~ 28 cm, 
respectively to derive Im[/x]. To produce the dots, we have assumed that complex Y v is the 
only source of CP-violation and have randomly produced Y v compatible with the data. The 
blue dashed line is the present bound on d e [18] and dot-dashed lines show the values of d e 
that can be probed in the future [19, 20] 
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